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Abstract 

An exact solution of a two dimensional RSOS model of wetting at a corrugated (periodic) 
wall is found using transfer matrix techniques. In contrast to mean- field analysis of 
the same problem the wetting transition remains second-order and occurs at a lower 
temperature than that of the planar system. Comparison with numerical studies and 
other analytical approaches is made. 

Introduction 

The study of d dimensional inhomogeneous fluids has traditionally focused on 
systems which retain some degree of translational invariance |l[|. Consider, for 
example, the adsorption of a fluid at a solid substrate modelled as a structureless 
wall which exerts a one-body potential on the fluid particles. If the interface is 
planar the equilibrium local number density depends on the normal distance from 
the wall (z, say) and the excess free energy is proportional to the area Ar If 
we suppose that the bulk fluid (infinitely far away from the wall) corresponds to 
a phase a (e. g. vapour) attention focuses on the possibility of adsorbing a layer 
of phase (3 (e. g. liquid) at the wall-a interface. At bulk two-phase coexistence 
one of two scenarios is possible. Above the wetting temperature T n the thickness 
of the adsorbed layer is infinite while below T n it is finite. The wetting phase 
transition occurs as T — > T~ and may be first- or second-order corresponding 
to the discontinuous or continuous divergence in the height of the a (3 interface 
(measured from the wall), respectively. 

For a non-planar or rough boundary the problem poses new difficulties due to 
an even greater loss of translational invariance. However, recently a number of 
independent approaches have been developed to meet such a challenge |2], [| |, ||] . 
Simple phenomenological arguments P, |3| indicate that the increase in contact 
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angle 9 [I] at the non-planar interface is controlled by the roughness parameter 



where A is the surface area of the rough wall. While this has been confirmed 
at low temperatures T <C T.,, using rigorous methods [[|, mean- field studies [Q] 
suggest a more complicated scenario in the vicinity of a planar second-order wetting 
transition temperature. Consider, for example, a corrugated wall described by a 
local height variable zw{x,y) = V2asinqx where a is the root mean square width 
of the corrugations and — is their period. Analytical analysis based on a Landau- 
type free energy functional || shows that the transition is driven first-order for 
a > a*, where a* is smaller than the bulk correlation length of the (3 phase. 
This first-order wetting transition occurs at a lower temperature than the second- 
order phase boundary (denoted T n ) for the planar system. For a < a* the phase 
transition remains second-order and occurs at the same temperature T n , i. e. there 
is no wetting temperature shift. 

It is natural to question the robustness of this result when fluctuation effects, 
occurring outside the mean-field approximation, are included. In d = 3 (which 
is the marginal dimension for short-range forces) one would expect the topology 
of the surface phase diagram to remain largely unchanged ||. However in d = 2 
fluctuation effects may be extremely strong and modifications are quite possible. 
Indeed, recent numerical work suggests that the important parameter is not r but 
Cs, the roughness exponent for the substrate. Clearly, for a corrugated wall (s — 
but Giugliarelli and Stella H have considered self-affine substrates for which the 
average height fluctuation for a sample of length L scales like gl ~ L^ s . Their 
numerical results strongly indicate that for (s < (, with £ = ~ the thermodynamic 
roughness exponent 0, the transition remains continuous. However for > £ it 
appears to be first-order. 

In the present article we describe details of an exact transfer matrix analysis of 
wetting at a corrugated wall which confirms the numerically based suggestion that 
the transition remains second-order. As we shall see the analysis is considerably 
more involved than the planar example and requires some care. Importantly we 
shall show that in contrast to mean-field expectations the wetting transition tem- 
perature is reduced by corrugation even if the order of the transition is unaffected. 

Consider a wall corrugated so that even sites are a lattice spacings higher than 
odd sites (see Fig. 1). Our approach evokes a restricted solid-on-solid (RSOS) ap- 
proximation to the semi-infinite Ising model (valid at low temperatures) which has 
proven to be very profitable for examining wetting transitions [0. Any overhangs 
of the interface are ignored (as are bubbles of a phase in the j3 phase and vice 
versa) so that a single valued variable measuring the height of the interface from 
the wall can be used. To finish our introduction we briefly recall the Hamiltonian 
for the planar system (a = 0) and suppose that the wall has N sites and contact 
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Interface 



U2n-1 



Site 2n-1 
Corrugated Wall 

Figure 1: Schematic illustration of a corrugated wall of amplitude a. A typical 
interfacial configuration (dashed line) is also shown. 



potential W > 0, 



N 



H[ u \ = Yl { J \ u n+i - u n \ - W5 Unfi } 



(2) 



n=l 



The variable u n measure the height of the interface at site n and J > 0, is the 
usual coupling parameter which suppresses bending fluctuations. Here periodic 
boundary conditions imply that un+i = u\ and with the RSOS constraint adjacent 
heights are restricted such that \u n +i — u n \ G {0, 1}. With this assumption we note 
that the Boltzmann weight is 



-J\u u |/T _ j(§ u , u+1 _|_ 5 u ' jU ^i) + 5 UtU > 



(3) 



where we have introduced the notation j = e J I T (with < j < 1) and set ks = 1. 
The model exhibits a second-order wetting transition when w = e w ^ T satisfies M 



w 



l + 2j 
1+J 



(4) 



which is the correct low temperature limit of the exact Ising model result 



RSOS model with a corrugated wall 

As adjacent interface configurations are restricted there are only two cases that 
need to be considered 

(i) a = 1. For this value of a the interface can access every site of the wall and 
the Hamiltonian is 

N 

H\[u] = Y {A U 2n - U 2n -l\ + J\u 2n - U 2n +l\ - W / (5 M2n _ lja __i + 5 U2nja )| (5) 
n=l 
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(ii) a > 1. For all a > 1 the interface cannot reach odd sites and still obey 
\u n+ i — u n \ G {0,1} (see Fig. 1). In other words there exists sites that are 
never accessed. While this may appear to be an artifact of the RSOS con- 
straint, in an unrestricted interface such configurations would cost enormous 
bending energy and so would be very rare. Consequently we expect the 
RSOS approximation to produce little error, particularly for larger values of 
a. The Hamiltonian reads 

TV 

#>lM = H { J \ u 2n ~ «2n-l| + A u 2n ~ U 2n +l \ ~ W6 U2n>a } (6) 
n=l 

It is in fact possible to account for both these scenarios with only one Hamil- 
tonian 

N 

H \ U \ = {j\ U 2n ~ U 2n -l\ + J\u 2n ~ U 2n +l \ ~ W*8 U2n _ lta ^ 1 - W5 U2nia j (7) 

n=l 

where here there are 2N sites along the wall, u 2n > a, u 2n -\ > a — 1 and periodic 
boundary conditions are still imposed (u 2 n+i — u i)- F° r a = 1 and a > 1 the 
effective variable W* takes the values W and 0, respectively. 
The partition function Z may be written as 

N 

/(«2n-l,«2n,M2n + l) I Q\ 



{u} n=l 



where / is defined to be 



f{u, u', v) = J\u - u'\ + J\v! -v\- -W* {8 u<a ^ + S v%a - X ) - W5 u , ta (9) 

To simplify (H) we follow 0] an d write 

Z= J2 5 P , Ul e~ f{u2N ~ 1,U2N ' Up) ■ .. e - f{ui ' U2 ' u:i) (10) 

{u},p 

Introducing a complete set of normalised eigenfunctions such that 

E$°1$ = «n.m (11) 
i 

with the over-bar denoting complex conjugation, allows Z to be rewritten 

Z= ^2 i(j^e~^ U2N ~ 1 ' U2N ' UpS> ■ ■ ■ e~f( Ul ' U2 ' U3 " l 'ip^ (12) 

i,p,{u} 

If the eigenfunctions are chosen to satisfy 

oo 

]T e^"'^^ = A (i) 4 i} (13) 



u=a— 1 
u'=a 
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then by using ( pT3|) repeatedly, (|T2| ) becomes 



Z 



iV 



l.p 



E(A« i| ) 
(A<»>)' 



A? 



'AW' 
A(°) 



A? 



+ ■ 



where A*- - 1 is the maximum eigenvalue. In the thermodynamic limit iV 
free energy density f reduces to 



(14) 



oo the 



f 



lim 

iV— »oo 



log( 



logZ 
2N 



(15) 



where we have used (p!4|). As usual the calculation of the free energy reduces to 
the problem of finding the largest eigenvalue of the matrix equation (|T3"D. 



Solution of the eigenvalue equation 

Consider (|T3| ) when i = 



(16) 



u=a— 1 
u'=a 



where the superscripts have been dropped for convenience. Using the Boltzmann 
factor (131) this reduces to 



E Wu 

(^a-l+lVa-l) 5 , 



u=a— 1 
u'=a 



(17) 



For values of u > a + 1 , ipv is beyond the influence of the wall and (|T7]) becomes 
fifiv-2 + 2jW-i + (1 + 2j 2 - \)ifj v + 2j^+i + f^ v+2 = (18) 
The linearity suggests the ansatz ip v = t v for some constant t. This gives 

t A + 2j~H 3 + (j- 2 + 2 - \r 2 )t 2 + 2j-H + 1 = (19) 
which in turn simplifies to 



t 2 + - V\)t + ll It 2 + + V\)t + i 



o 



(20) 
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Denoting the roots t\, t% and £3, £4 for each bracket we see that 



and 
and 



ti + t 2 = -r 1 (i-v / A) 
h + u = -r x (i + VX) 



with (|i~8| ) having a general solution 



^ = Atfl + B t v 2 + A x tl + B x t\ 



(21) 



(22) 



for some constants Ao, A\, Bq and B\. However we expect the eigenfunction ip. v to 
be finite as v — > 00 and so if the are real we can choose without loss of generality, 
see ([H]), |^3| < 1 and \t2\1 \t^\ > 1, which implies that B and Bi must vanish 
in (p2|). Notice that this is not true if the roots are complex and then ti—t^ — e %e 
for some 8. By considering the discriminants of fl2"0|) it is straight forward to show 

(23) 



tl 


= -e- p ; y/X 


= 1 — 2j cosh 


; v^A < 1 - 2 j 


h 


= e id ; V\ 


= l + 2j cos 


; 1 - 2j < y/X < 1 + 2 j 


ti 


= e"^ ; VA 


= 1 + 2j cosh /i 


; > 1 + 2j 


,P 


and 6, and 








t 3 = -e-' il ; 


\/A = 2j cosh 


- 1 ; v^A > 2j - 1 




ti = e idl ■ 


vA = 2j cos — 


1 ; VA < 2j - 1 



for some /ii and #i. The eigenfunction ^ is 

tp v = Aotl + A,t v 3 



(24) 



(25) 



for v > a + 1. In the thermodynamic limit only the maximum eigenvalue con- 
tributes to the free energy and so for low temperatures (for which the eigenfunc- 
tion must correspond to a bound state) we have v^A > 1 + 2j > 2j — 1 and ( [25]) is 
valid for 



h 



- 2 \3 



\Vx 



i( v ^_l)2_ 4 | 



-Mi 



r-2 



(26) 
(27) 



(VA + l) 2 -4 

and ( P4[ ) this will happen 

(28) 



To deal with the boundary conditions, i. e. (|ITD for values of v < a + 2, we extend 
the solution (|25|) to the range w > a + 1 and impose 



Wetting will occur when any ti turns complex. From 
first for ti when v^A decreases to 

VX = 1 + 2 j 



for all v < a — 1 



(29) 



6 



due to the RSOS constraint. Then ( [P7|) for v — a — 1, a, a + 1 and a + 2 gives a 
set of equations for ipa-i^aj A) an d A- These are best written in matrix form 

/ j 2 ww* - A jw^fw* j 2 t 2 w^fw* j 2 t 2 wy/w* \ 

jwVw* j 2 + w - A jt\ (1 + jti + to) + i*3 + 

i^v 7 ^ + jti[(w;-l)iti-2]-i 2 jt 3 [(^-l)ji 3 -2]-i 2 

V j 2 -fh -j 2 t 3 

( V>a-1 \ 

A 

V M J 

To have non-trivial solutions the matrix, M. say, must have zero determinant 



(30) 



det.M = (31) 

which is a (complicated) equation for A in terms of t\ and t$. From here one 
can of course proceed analytically but the task is laborious and ideally suited for 
computer algebra packages. Using MATHEMATICA A can be obtained exactly 
and then used to solve (|30|), remembering the normalization condition Yl v l^vl 2 — 1- 
The resulting expressions for if) a -ii i^a, A and A\ are very long and not particularly 
illuminating. 

However, exactly at wetting holds and it is easy to demonstrate from (|26|) 
and that 

ti = 1 (32) 



* 3 = -(r 1 +i)+r 1 v 1+2 J ( 33 ) 

and defining 



q = yjl + 2j ^Kq<VS (34) 
equation (|3l|) can be shown to become 

(q 2 - l) 2 w*w + 2q(l + q - q 2 + q 3 )w = 4g 4 (35) 

which gives the wetting temperature. 

As an aside, we point out that our method of solution is quite versatile and can 
be very simply adapted to study the situation considered by Nechaev and Zhang 



11]. These authors have a periodic potential at a planar wall — W taking the 



values Wo and W\ at odd and even sites respectively. Their Hamiltonian is 

N 

H \ U \ = { J \ U 2n ~ U 2n -l\ + A U 2n ~ «2n+l| ~ ^0^-1,0 ~ Wl^.o} (36) 

n=l 

which is very similar to (|7|). Indeed, away from the wall still holds and only 
the elements in the matrix Ai change. Wetting again occurs when vA = 1 + 2j 
and using the same techniques as before the wetting temperature is found to obey 

4g 4 + (1 + q 2 )(q 2 - 2q - l)w w 1 + 2q 3 (l - q)(w + w x ) = (37) 
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where wq = e w °^ T and w± = e Wl ^ T . With a little effort this can be shown to be 
equivalent to equation (14) in which provides a useful check for the validity of 
our method. 



Comparison of wetting temperatures 

We consider the two possible wall configurations separately 



(i) a = 1. In this case w* = w and (|35|) is a quadratic equation in w. Hence we 
arrive at 



-q 1 -qf + ql-qi + q 1 Jl + 2q 1 + 3q 2 1 -5qt-2ql + 5qf 
Mil) = (g 2 _ X y ( 38 ) 

for the wetting phase boundary w\ = e w ' Tl . For fixed values of W and J 
the wetting temperature T\ is reduced from its value T n in the planar case, 
which from (01) satisfies 



2q 2 

W ^=lh (39) 

To see this first assume 

T x > T v (40) 
Then since w = e w ^ T and j = e~ J//T we have 

w n > Wi (41) 
qi > q* (42) 

as 3i > j-n- Using (|38| ) and fl39]), equation (|4|) implies the inequality 

(1 - gi ) gi (l + qi ) 2 ql + (1 + q x - q\ - qf + q$)q$ > q\ (43) 

or (as qf > g*) 

(1 - qi ) qi (l + qi f + gi(l - gi) 2 (l + qi)ql > (44) 

This then gives 

1 + qi + (1 - < (45) 
but q 2 < 3 from fl3"4] ) and so 

l + gi + 3(l- gi ) <0 (46) 



indicating that q\ > 2. However, this is a contradiction — from (|34] ) we have 
q < a/3 and hence the assumption (|40| ) is false. Thus Ti is lower than T % . 
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(ii) a > 1. Now W* = implying w* — 1 and (p5|) gives 



^ = i + 2 92 -ti + 3,i (47) 

where we write u>2 = e^ 72 denning the wetting temperature T 2 for a > 1. 
The shift of the wetting phase boundary can be seen following a similar 
method to that given above. 

First assume 

T 2 > 71 (48) 

which then implies, 

wi(qi) - w 2 (q 2 ) > (49) 
q 2 > qi > 1 (50) 



and recall from ([34]), 1 < q 2 < V3. It is straight forward to see, via 
that w' 2 (q) > for all 1 < q < and u>2(l) = 1- Hence (|50"1) gives 
< w 2 (qi) < w 2 (q 2 ). Then, defining w(q) = W\(q) — w 2 (q), fl49^) becomes 

w(q 1 ) > Wi(qi) - w 2 (q 2 ) > for 1 < q x < (51) 

However one can also show from (|38D and (0) that w'(9) < for 1 < q < \/3, 
and tf(l) = 0. Hence w(q) < for this range of q which contradicts (]5l"l) 
based on the assumption 



In summary we have proven that 

T 2 < T x < T w (52) 

so that the larger the magnitude of corrugations of the wall the further reduced 
is the wetting temperature. Inequality (|52"D is simply explained using entropy 
arguments in the conclusion. 

Finally we turn to a calculation of the critical exponents characterising the 
transition in the non-planar system. First consider the mean interface position 
above the wall defined by 



E{n} ue 



-HU 



^ } „_ g , M i (53) 



Em 



In the thermodynamic limit this can be written (for any odd numbered site) 

{uUd ~ EJ^P (54) 

where we have expressed the wave-function normalization explicitly. The sums in 
([54]) are simple geometric progressions and are easily evaluated. Expanding the 
result around the critical point /i = shows the leading-order contribution to be 

i w — (55) 
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where \x ~ Tj — T and so the adsorption critical exponent /? s = 1 identical to the 
planar case JF]. Clearly the phase transition remains continuous and falls in the 
same universality class as the planar system. 

Similar remarks apply to the singular contribution to the excess free energy 
and transverse correlation length £y. For example, for the free energy density we 
find 

f = -logVA 

-log(l + 2j)- r ^- (56) 

for small /i, identifying the specific heat exponent a s = 0. Again, the exponent 
describing the divergence of the correlation length v\\ = 2 as expected 0. 



Conclusion 

We have presented an exact analytic solution to the problem of wetting in two 
dimensions with a corrugated wall. The corrugations can be of arbitrary height 
without the transition losing its second-order character. The critical exponents 
governing thermodynamic quantities as the phase boundary is approached are un- 
changed from the planar case although the wetting temperature itself is reduced. 
These results contrast sharply with mean-field theory [[| but support numerical 
studies in d — 2 with self-affine (rough) walls []5). In fact both the d = 2 and 
mean-field result suggest that by including thermal fluctuations second-order wet- 
ting transitions become roughness-induced first-order if the width of the wall (as 
measured by a roughness exponent (s for d < 3 and length scale a for d > 3) is 
larger than the analogous quantity for the fluid (a(3) interface. 

Our explicit calculations illustrate that the transition temperature is reduced 
by wall corrugation and to complete our article we make some comment on the 
generality of the results. Below the upper critical dimension wetting transitions 
are bought about by entropy winning its competition with internal energy in the 
free energy of the interface. For a corrugated wall the internal energy of a bound 
state is less negative than the planar case because either there are fewer points 
of contact with the substrate if the a(3 interface is flat or an increased bending 
energy if it leaves the plane to reach the contact points. Consequently entropy, 
which is hardly affected by the corrugation, needs to overcome an internal energy 
of reduced magnitude implying that the wetting temperature is decreased. 

With this general argument we expect that the RSOS result described here 
captures all the essential physics of the full Ising model calculation (at least at low 
temperatures). While an Ising model solution would be very welcome the transfer 
matrix analysis would be much more complex and may be prohibitively difficult. 

We acknowledge financial support from the Engineering and Physical Science 
Research Council, United Kingdom. 
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